Numerical Test of Disk Trial Wave function for Half-Filled Landau Level 



S.-R. Eric Yang 
Department of Physics, Korea University, Seoul 136-701, Korea 
Asia Pacific Center for Theoretical Physics, Seoul, Korea 

Min-Chul Cha 

Department of Physics, Hanyang University, Ansan 425-791, Korea 

Jung Hoon Han 

Department of Physics, U.C. Berkeley, Berkeley, CA 94720 USA 

The analyticity of the lowest Landau level wave functions and the relation between filling factor 
and the total angular momentum severely limits the possible forms of trial wave functions of a disk 
of electrons subject to a strong perpendicular magnetic field. For N, the number of electrons, up to 
12 we have tested these disk trial wave functions for the half filled Landau level using Monte Carlo 
and exact diagonalization methods. The agreement between the results for the occupation numbers 
and ground state energies obtained from these two methods is excellent. We have also compared 
the profile of the occupation number near the edge with that obtained from a field-theoretical 
method. The results give qualitatively identical edge profiles. Experimental consequences are briefly 
discussed. 



I. INTRODUCTION 



There has been a great interest in two-dimensional electron liquid in a strong raastietic field at half-filling. Its bulk 
properties have been investigated and several novel features have been revealedu aaa. Electrons form a new exotic 
metal in wn ifbiAbpy move effectively in zero magnetic field despite the presence of an external fieldo. A dipole-like 
new particleoQU'El with, fcrmionic statistics has been introduced. This new state of matter forms a Fermi liquid with 
a finite compressibility!! p. n I— . 

The edge properties of this exotic liquid, are of great current interest, both experimentallytl and theoreticallycTo. 
Recently a two-mode model for the edgalj has beeruproposed based on a theory of the electron liquid at v = 1/2 
consisting of charged Bose and neutral Fermi liquidsa. In this theory the charge mode propagates while the neutral 
mode is dispersion-less. A theory based on the topology of Abelian fractional quantum Hall states also leads to similar 
results, but in this approach two non-propagating modes of a purely topological origin existEl. 

A disk of electrons consists of a bulk and an edge, and it is desirable to treat them in the same theoretical framework. 
The edge-and bulk properties have been also investigated using a trial wave function for an electron disk at half- 
fillingtL3ii-!l. The trial wave function was constructed using analyticity of the lowest Landau level wave functions, and 
the fact that the filling factor (v) is related to the total angular nacpacntum (M z ) through v = N(N — 1)/2M Z . These 
conditions severely limit the possible forms of the wave fimctionOlij. It was tested that the wave function overlap of 
this state with the exact one is essentially one for N = 50. The computed edge occupation number showed, in sharp 
contrast to incompressible Lpjghlin states, that a tail exists in the occupation number beyond the radius of a_droplet 
of uniform filling factor l/2t3. The properties of zeroes of the wave functions have been also investigateefcj. The 
proposed trial wave function can be used to understand the dynamic structure factor, which is a physically measurable 
quantity. This was done by constructing trial wave functions for some of the low energy excitations of a diskLj. One of 
them is the center of mass (CM) motion, which can be interpreted as an edge mode. Other excitations are identified 
and their wave functions are also constructed. They are not edge modes, but are excitations of the entire droplet. 
They are new excitations of a droplet at v = 1/2 which do not exist in incompressible states. The dynamic structure 
factor shows that the disk can support several low energy excitations, unlike incompressible quantum Hall states. 

In this paper we perform an extensive numerical test of the proposed-trial wave functions for an electron disk at 
half- filling. We find that the previously proposed trial wave functionald'E^I work quite well up to N = 8. We find that 
for bigger values of N almost degenerate states with lower energies exist. Using these trial wave functions we have 
computed several physical properties by performing Monte Carlo simulations, and have tested the obtained results 
against exact diagonalization results for up to N = 12. The agreement between the results for the occupation numbers 
and ground state energies obtained from these two methods is excellent Wp have also compared the profile of the 
occupation number near the edge with that of a field theoretical approachlHJuj. They give qualitatively identical edge 
profiles. 
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This paper is organized as follows. In Sec. II we explain the composite Fermion wave function, which is a central 
concept in our trial wave function approach. Next, in Sec. Ill, we give the physical idea behind the construction of 
trial wave functions at the principal filling factors using composite Fermion wave functions. In Sec. IV we perform 
Monte Carlo simulations using the proposed trial wave functions and compared the extracted physical properties 
with those of exact diagonalization method. A field theoretical calculation of the edge occupation number is carried 
out in Sec. V and the obtained result is compared with those of the exact diagonalization and Monte Carlo results. 



Discussions are given in Sec. VI. In this paper length and energy are measured in units of > 
magnetic length. 



and e /e£, where £ is the 



II. COMPOSITE FERMION WAVE FUNCTIONS 

In this section we explain composite fermion wave function, which is a central concept in our trial wave function 
approach. This concept was introduced by Jain and Kamilalij. Before we discuss our trial wave function approach it 
is necessary to understand this concept first. So we give a short summary of their result. 

A. Physical Model of a Disk 

Let us first set up our model of a disk. We consider a system of 2D electrons confined by a parabolic potential, 
V(r) — m e 51 2 r 2 /2. The advantage of this model is that some exactly known results can be utilized in constructing the 
collective modes. Such a system will form a uniform density electron disk for sufficiently large N. The symmetric gauge 
single-particle eigenstates 4> n ,m(z) are conveniently classified by a Landau level index n and an angular momentum 
index m — — n, — n + 1, ... (z — x + iy). In the strong field limit the single particle orbitals in the n = level 
have energies e m — Huj c /2 + j(m + 1), where 7 = m e £l 2 £ 2 . The system is invariant under spatial rotations about 
an axis perpendicular to the 2D plane and passing through the center of the dot. It follows that the total angular 
momentum M z is a good quantum number. Eigenenergies may be expressed as a sum of interaction and single-particle 
contributions 



Ei(N, M z ) = Ui(N, M z ) + 7(JV + M z ). 



(1) 



Here i labels a state within a M z subspace, and Ui(N, M z ) oc e 2 /e£ is determined by exactly diagonalizing the electron- 
electron interaction term in the Hamiltonian within this subspace. We will present in this paper a different approach 
than some of previous works on fractional quantum Hall states disk geometrytSEH. 



B. Composite Fermion Wave functions 



Recently it has been argued that the-jpain physics of strong electron correlation in quantum Hall systems is that 
the electrons carry two Jastrow factorsES. 



iI>[z] = Pl* Y[(zj-Zk) 3 

j<k 



0mmi(zi) <f>mmi{ZN) 



(2) 



where P^3 stands for projection to the lowest Landau level (LLL). It is useful to rewrite this as 

(z N )J N 

i/>[z] = P L 3 ' 



r (zn)Jn 



(3) 



In physical terms, the Jastrow factor J; = Jlfc^j( z « — z fe) implies that each particle carries with it a correlation hole 
with respect to all the other particles. ■— ■ 

The projection P L 3 is tantamount to replacing the anti-holomorphic friable z with 2(d/dzp3. It is not straight- 
forward to project this wave function to the LLL. Jain and Kamila (JK)tj have proposed the following wave functions 
of interacting fermions in the LLL for a given total angular momentum M z : 



(zn) 

r )n N m N {,Zi) 1]n N m N (zN) 



(4) 



2 



d™Ji 



(5) 



Vnm{z) is not a genuine single-particle state as its value depends on the position of all the other particles. Nonetheless 
to stress the similarity between Eq.(^|) and Eq.(|J)JK regard n and m as the LL index and the angular momentum of 
the composite fermions (CFs), respectively. The composite fermion Landau level (CFLL) index is just a convenient 
nomenclature for some abstract mathematical property of rj nm (z). It should be noted that CFLL index n is not a real 
Landau level index, and although higher CFLL indices appear in this wave function it lies entirely in the lowest electron 
Landau level. On the other hand, it can be shown that the total angular momentum of ^cf is M z = N(N — 1) + M*, 
where M* = J^. m; is the sum over angular momentum of CF levels (N(N — 1) stems from ./V Jastrow factors of Jj). 
JK choose the ground state as the compact state [Nq,Ni, .., JV/_i], where the number of CFs in the i-th CFLL is 



Ni. This state is defined to satisfy: [1] N > Ni > N 2 , [2] ELo^ = N > I 3 ] For S iven N i> a11 the CF levels 
m = —i, —i + 1, .., Ni — i — 1 are assumed to be occupied. For example, for N = 5 the CF occupation numbers are 
given by (•, •, •, •, •, o, o, ...) in the i-th CFLL [An occupied (unoccupied) CF state is represented by a full (empty) 
circle. The angular momentum of a CF state increases from left to right]. The ultimate justification of choosing 
these wave functions comes from the comparison with exact diagonalization results (See Sec. IV.B.l Wave function 
overlaps) . 

We will drop the universal factor FL e~' Zi / 4 from now on. When we write out the wave function of [No, AT 1; .., iVjj i] 

it has the form 



Ji 

J-2 



Z\J\ 
Z 2 J 2 
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Jn z^Jn 



^1 diJi 
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(6) 



Successively higher powers of di acting on Ji appear as we move to th«ight for a given row. The matrix elements are of 

the form z" + " l <9" Ji. One can easily verify that the 1/3 Laughlin statecj is given, in this notation, by [No, AT 1; .., iVjj 1] = 

[7V,0,...,0]: 



ZxJx 



,JV-1 



■h 



Jn zmJn z^ 1 Jn 



(7) 



It can also be verified that the wave function for the filled electronic Landau level at v = 1 is given by 
[JVo.-tfi,..,^-!] = [1,1,..,!]: 



Ji d\J\ 
Jn 9nJn 



1 Jn 



(8) 



JK do not give specific prescription for other general filling factors. 



III. CONSTRUCTION OF TRIAL WAVE FUNCTIONS FOR PRINCIPLE FILLING FACTORS 



Now that we have explained the composite fermion wave functions we can proceed to construct trial wave functions 
at other principal filling factors. 



A. Uniform density droplet 



The formula v = N(N — 1)/2M Z , which relates-the filling factor to the total angular momentum can be shown to 
be true for Laughlin states using plasma analogycJ. For other states described by composite fermion wave functions 
numerical results support the validity of this formula. We have tested numerically the wave functions [4, 3, 3, 1] and 
[6,4,3,2,1,1], which we believe describe the v — 1/2 state at N = 11 and N — 17, respectively. We have verified 
for these states that the average electron density n(r) inside dot is approximately (l/27r£ 2 )^ (see Fig. [j]). We have 
also tested other states [6, 1, 1] and [9, 1, 1] for N = 11 and N = 17 and found similar results. We expect from 
the size dependence that the agreement should become better for N >> 17 (17 is the biggest size we could handle 
numerically). Even for the Laughlin states numerical results show that the correct average value is approached slowly 
as N increasestil. 
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B. First method 



Let us first explain our trial wave function approach using analyticity of the lowest Landau level wave functions 
and the fact that the filling factor is related to the total angular mommtum. 

As a trial wave function for the state at v = 1/2 Yang and HanOO have proposed the state [N/2, 2, 1, 1] 
([(N — l)/2, 1, 1]) for even (odd) N. It can be verified that for each value of N = 3, 4, 5, 6, 7 this state is the only 
possible compact state satisfying the condition v = N(N — 1)/2M Z . The occupation of CFs in these states are listed 
in Table § The state [f , 2, 1, 1] has the form 
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This state has the following form in terms of antiholomorphic variable z 



Ji 



— -i 



Jl ZlJl 



— -1 



Jl 



Jn z^Jn ^jv-i Jn z^Jn z N Jn 

Han and YangE3 have also proposed trial wave functions at other principal filling factors. 



(10) 



C. Second method 



The proposed wave functions can be also constructed using a more physical argument based on counting flux 
numbers and their relation to the total angular momentum. 

Consider the v = 1/3 state [N], which has 3N flux quanta. The angular momentum of this state is 3N(N — l)/2. 
If two flux quanta are removed from this state the number of flux quanta per particle will be 1/z/ = ,7 2 . Then, in 
the large TV limit, the total angular momentum will be given by 



N 



M z = ^N(N - l)/2 = 3N{N - l)/2 - (N - 1) « 3N(N - l)/2 - N. 



(11) 



^From this we see that reducing two flux quanta amounts to reducing the total angular momentum by N. The 
resulting state is [N — 1, 1] since the angular momentum of this state is N(N - 1) + M* = 3N(N — l)/2 — N. 
Here we have used the fact that the total angular momentum of composite fermions M* in the state [N — 1,1] is 
— 1 + X)ilo 2 m i = ~ 1 + {N ~ 1)(-W — 2)/2. If two flux quanta are removed additionally, the angular momentum is 
again reduced by N and the resulting state is [N — 2, 1, 1]. If two flux are removed q times successively the resulting 
state will be [N — q, 1, 1, 1], where 1 appears q times (However, this state is not unique. For other possibilities, sec 
below). The total angular momentum of this state is 3N(N - l)/2 - Nq. The half filled state has M z = N(N - 1) 
so q = (N - l)/2. For the v = 2/3 state M z = 3N(N - l)/4 so q = 3(N - l)/2. For v = n/(2n ± 1) we find 



_ (N-l) 



(3-^). 



D. Trial wave functions for A*' > 9 



As the number of electrons increases, other choices of occupations are allowed than proposed previ- 
ously. For v = 1/2, unlike the case for N < 7, there are two possible compact states at N = 8 
with the correct value of M z = N(N - 1) + M*\ The first state is [4,2,1,1] with occupied CF states 
(n,m) = (0,0), (0,1), (0,2), (0,3), (1,-1), (1,0), (2, -2), (3, -3), and the second state is [3,3,2] with (n,m) = 
(0,0), (0,1), (0,2), (1,-1), (1,0), (1,1), (2, -2), (2,-1) (in both cases J2 t m i = M * = °)- Note that the second state 
does not have the proposed form [N/2, 2, 1, 1], but it also turns out its energy is almost degenerate with that of 
[4, 2, 1, 1] (the energies are 6.025 vs. 6.095). For v = 1/2, the possible choices of [No, N\, ■ ■ ■] are given in the Table | 
for several values of N. Possible choices [Nq, N\, ■ ■ ■] at other principal values of v are given in Table 0. The Monte 
Carlo energies of these states at v = 1/2 are given in Table |[ We choose the state with the lowest energy as the 
groundstate. 



4 



IV. MONTE CARLO AND EXACT DIAGONALIZATION RESULTS 



Ground state properties can be calculated from the proposed trial wave functions. Here we concentrate on occupa- 
tion numbers and ground state energies. 



A. Monte Carlo calculation of occupation numbers 

For Laughlin states the edge occupation number n(k) ~ (fe — fc e d)' ' u ~ In disk geometry single particle 

angular momentum m is related to the wave vector through, k = 2iim/ R m , where R m = y/2(m + 1). The value of m 
corresponding to the location of the edge is m ec i = 2N — lE3 for v = 1/2. The corresponding k is defined as k ee [. The 
exponent, 1/v, in this expression is intimately related to the short range correlations of the Laughlin liquid. Therefore 
one can expect to gain information about the the short range correlations from the edge occupation numbers. This 
is the motivation for investigating occupation numbers. 

Now we evaluate the one-particle density matrix n{z, z'), which is defined as 

n(z,z') = N d 2 z 2 ■ ■ -d 2 z N ip*(z 7 z 2l ■ ■ ■ , z N )ip(z' , z 2 , • • • ,z N )/Q N , (12) 



where the normalization integral is Qn = (ip\ip). We are especially interested in the angular dependence of the 
one-particle density matrix, n(r 1 re l6 ) 1 where r = \z\. Using n(r,re l9 ) — J d 2 zi(l/2Tr\zi\)S(\zi\ — r)n(zi, zie lS ), we 
have 

n(r, re ie ) =J^f d ' ^ ' ' ' ^ z n S ^\ ~^ 4>*(zi, ' • ■ , *i, ' ' • , *atM*i, • • • , z^\ z N )/Q N , (13) 

i ' 1 

It is convenient to perform the integration by Monte Carlo sampling with the probability distribution 

P{z\,z 2 , ■ ■ ■ ,z N ) = \ip(z 1 ,z 2 , ■ ■ ■ ,z N )\ 2 /(ip\ip), so that 

n(r,re ) = ) d z\ ■ ■ ■ d z N — — — , — P(z 1 , z 2 , ■ ■ ■ , z N ). (14) 

~^ J 2tt\z 1 \ ip(zi,z 2 , ■■■ ,Zi,--- ,z N ) 

The angular-momentum occupation number is determinecS by 



1 ^ e^ m n(r,re^) 

where m e ft is the highest power of Zj, and 



r , Q meff + l ' 1 } 



is 



u (r)ia = ^- r — e -r 2 /\ (16) 
l ^ K n 2nm\ 2 m V ; 

The interaction energy can also be estimated from these wave functions. The interacting part of the the Hamiltonian 

(17) 

Monte Carlo sampling has been carried out through Metropolis algorithm by changing the position of electrons, 
z.i — > Zj + Szi, where typically \8zi\ < 0.2 — 0.4. The wave functions vary rapidly with respect to the positions of 
electrons so that after one Monte Carlo sweep IV'Inew/lV'lold cnan g es U P to ~ 10 10 or 10~ 10 . It is remarkable that, in 
spite of this rapid variation of the wave functions, Monte Carlo method quite nicely gives accurate expectation values. 



5 



B. Numerical Results 



1. Wave function overlaps 

The justification of the proposed wave functions at v = 1/2 comes from small-size exact diagonalization. For 
N = 5 apd for the Coulomb interaction the overlap of this ground state wave function with the exact wave function 
is 0.9989O. (The accuracy of the proposed trial wave function is comparable to that of the Laughlin state: for four 
particles with Coulombic repulsions the overlap is 0.979 for the v — 1/3 Laughlin stateH3). For N — 5 trial wave 
functions of the first, second, and third excited states with the change in the total angular momentum 5M Z = 1 can 
be constructed from the ground state wavefunction. The overlaps with the exact ones are 0.9990, 0.9663, and 0.9098, 
respectively^. 



2. Occupation numbers and ground state energies 

For N bigger than 5 the proposed wave function of Eq. (|]) can be tested against the exact diagonalization results 
by comparing the occupation numbers. 

Fig. |2| displays occupation numbers for v = 1/3 disk. The results obtained from the Laughlin wave function and 
those obtained from exact diagonalizations are shown. The agreement between the two are excellent near the edge. 
This implies that the Laughlin wave function captures the short range correlations remarkably well. Note that near the 
center of the disk a small discrepancy exists between the results of Laughlin wave function and exact diagonalization. 

Fig. H displays occupation numbers for v = 1/2 disk, calculated using Eq.(^). Again, as in the case of the Laughlin 
wave function, the results obtained from the proposed trial wave function and those of the exact diagonalization are in 
excellent agreement. For N = 8 the state [3, 3, 2] has the correct angular momentum consistent with v = 1/2, but its 
occupation numbers show a significant deviation from the exact ones even near the edge. In contrast the state [4, 2, 1, 1], 
given by Eq.(^), agrees with the exact results. For N — 9 the energies for two possible wavefunctions [5, 1, 1, 1, 1] and 
[3, 3, 3] are nearly the same, deviating from the exact diagonalization result by about 1%. We find neither [5, 1, 1, 1, 1] 
nor [3, 3, 3] give accurate occupation numbers. In this case even the exact diagonalization procedure converges very 
slowly. Probably an elaborated linear combination of wave functions could lower the energy. For N — 10, 11, 12, 
we find, among possible choices of wave functions, that the results of the lowest energy states agree with the exact 
diagonalization result (see Table | and Fig. ||). Note that exact diagonalization results are not available for N > 13. 



3. Zeroes of wavefunction 

For TV = 12 the trial wavefunction [5, 3, 2, 1, 1] has the lowest Monte Carlo energy. We have investigated zeroes of 
this trial wavefunction (see Fig.^). We find that in [5,3,2,1,1] zeroes are closer to the electrons in comparison to 
[6, 2, 1, 1, 1, 1]. This is consistent with the fact that [5, 3, 2, 1, 1] has lower energy than [6, 2, 1, 1, 1, 1]. Properties of 
zeroes of other trial wavefunctions are discussed at length in Han and YangEJ. 



V. FIELD THEORETICAL CALCULATION OF OCCUPATION NUMBERS 

Sinali-size exact diagonalization calculation at v = 1/2 indicates that n(k) near the edge depends linearly on k for 
k < fc e J13. We investigate using a field theoretical approach whether this feature survives in the thermodynamicJimit. 
The hydrodynamic action for the edge mode of a half-filled droplet has been written down by Lee and Wer£3, and 
independently by Lopez and FradkinEH. Ignoring the coupling to the bulk degrees of freedom, the action is given by 

Sed3e = hjt x {t~ Vc t) dxdt ~ h I £ \Tt + Vn te) dxdt ' (18) 

where d x <fi and d x r\ are the displacement fields of the charge mode and the neutral mode. Following Lopez and 
Fradkin's approach one would have to add another neutral mode in the action. The results for the occupation number 
remain unaltered, and we choose to present our result with the two-mode model for simplicity. The boson fields <p{x) 
and rj(x) satisfy the following commutation relation 

= wiS ( x ~ V)'> [^ ,7? ( y )] = -2TriS(x - y). (19) 
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The electron operators are constructed ast2lliil 

t/j e (x) ~ e -2#(x)±<7 7 (x) > (20) 

One can verify that the operators satisfy the necessary anticommutation relations. 

The occupation number profile follows from the equal-time single-particle Green's function, G e (x) — (ipl(x)tp e (0)) , 
Since the dynamics of <j> and r] are independent, one obtains 

G e (x) ~ exp[G*(x) + G v (x)] (21) 

where 

' 1— exp[— 2tt(1— ixA c ) /A C L] " 



G<f,(x) = -2 In 



l-exp[-27r/A c L] 



G T) (x) = -lnp- CXp[ - 27r ; i+ ^" ) / A ^n (22) 
vK ' V l-exp[-27r/A„L] J V ' 

are equal-time Green's functions of <fi and to A finite circumference L of the droplet has been introduced, as well as 
some high momentum cut-off A c /A ra for charged/neutral modes. 

A Fourier transform of Eq. (|2l| ) should give the momentum occupation number. But it is first necessary to define 
the position of the edge. It is natural to define the position of the edge as m ec i from the the following simple model 
for the edge: Consider a disk of electrons at half filling with uniform occupation numbers given by n(m) — 1/2(0) for 
to < m e d{else). The sum over occupied n(m) is N so m e d = 2N — 1. A uniform disk at v = 1/3 has the occupation 
numbers given by n(m) = 1/3(0) for to < 3N — l(m > 3A*" — 1), implying that the edge is at m e d = — 1. We 
have used these values of m e d as the positions of the edges. In real systems effects neglected in this simple model 
will change the shape of occupation number around the edge. At half filling accurate numerical results demonstrate 
that the occupation number is smeared out around m e d — 2N — 1, just like at the edge of an ordinary metal. As 
Yang and Hanlla pointed out this definition of the location of the edge gives a good scaling curve (data collapse) for 
the edge profile for different values of N (other definitions will fail). In the v = 1/3 incompressible (non-metallic) 
case, the shape of occupation number around the edge is changed differently: it can be deduced from the Laughlin 
wavefunction that occupation numbers beyond m e d = 3N — 1 are zerccil. So in the incompressible (non-metallic) state 
of v = 1/3 the occupation number is zero outside the edge. We find that the momentum occupation number is 

C e -27rm/K„^ TO > 

n(m)/n(m = 0) = < „ , , ' . „ , (23) 
K " y ' \ e 2WKc{ 1 _ m ( 1 _ e -2*-( Kc +«„ ))} ; m < o. v ' 

Here K c , n is defined by A c ^ n L. The boundary of a compact v = 1/2 droplet coincides with to = (m is measured from 
the value corresponding to the location of the edge m e( j). We find that the occupation number falls off linearly with m 
for the interior of a compact droplet, and exponentially for the outside. Fig. [| shows momentum occupation numbers 
for eleven and twelve electrons at half-filling obtained from exact diagonalization as a function of x = (k — k e d)R e d 
(Red is the radius of the disk corresponding to m e d)- Monte Carlo results are shown for N = 13. These results show 
that electron occupation persists in the tail region, x > 0. They also show that near x < five consecutive data 
points fall, to a good approximation, on a straight line, and is well approximated by a linear behavior A + Bx for 
suitable choices of A and B, in agreement with our calculation, Eq. ([23). Based on the field-theoretical calculations 
presented above, we expect that this "occupation tail" will persist in the thermodynamic limit. 



VI. DISCUSSIONS 



Let us first give a short summary of our results. The analyticity of the lowest Landau level wave functions and the 
relation between filling factor and the total angular momentum severely limits the possible forms of ground state trial 
wave functions. The computed occupation numbers of these wave functions are in excellent agreement with exact 
diagonalization results. This was verified for N < 12 at v = 1/2. 

The proposed trial wave function at v = 1/2 can be tested against experiment. Low energy excited states can be 
constructed by following the prescription given in Yang and LyueEj. First the center of mass motion is constructed by 
multiplying the groundstate wave function with center of mass coordinate Z = Xi/N. Then the Gram-Schmidt 

process can be utilized to construct an excited state that is orthogonal to this state and the groundstate. Higher 
excited states can be constructed repeating the Gram-Schmidt process. Investigating small size systems indicate that 
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several low energy peaks are present in the dynamic structure factor of a disk of electronsE-l This feature is different 
from,that of the v = 1 incompressible state, which has only one low energy peak. Far-infrared radiation of a quantum 
dot)2a may be used to investigate these peaks. 

A naive incorporation of the linear dependence of the occupation numbers near the edge into chiral Luttinger theory 
suggests a value of tunneling exponent equal to 1. This value apparently agrees with experiment, but Fermi statistics 
will not be recovered in this naive approach. We believe that the tail in the occupation number beyond the edge must 
be incorporated to get the correct statistics. As far as the value of the tunneling exponent is concerned it seems that 
the important region is the linear region. 

For other principal filling factors beside v = 1/2 the edge profiles and dynamic structure factors can be investigated 
along the similar line developed in this paper. It remains a theoretical challenge to test whether the proposed trial 
wave function is equivalent to some of the field theoretical approaches developed recently. Test of the consistency of 
these different approaches will deepen our understanding the exotic half-filled Landau level. 
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TABLE I. Possible choices of [N ,N U - 


•] at v = 1/2 and energy evaluated. 




are the energies for each 


wave function. 


Eed are the energies obtained by the exact diagonalization method. Numbers in 


the parenthesis represent 


Monte Carlo error bars. 
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TABLE II. Possible choices of [No, JVi, • ■ ■] at general filling factors. 



N t/ = 2/5 v = 3/7 v = 3/5 v = 2/3 

3 

4 - [3,1] [2,1,1] 

5 [4,1] - - [2,1,1,1] 

6 - [4,2] [2,2,1,1] 

7 - [5,1,1] [3,1,...] 

8 [6,2] - - [2,2,1,...] 

9 [7,1,1] [6,2,1][5,4] [3,2,1,...][2,2,2,2,1] [3,1,...] 

10 - [7,1,1,1] [5,5] [4,1,...][2,2,2,2,2] 
11 

12 [9,2,1] [8,2,1,1] [4,2,1,...] [3,2,1,...] 

13 [10,1,1,1] [9,1,...] [5,1,...] [4,1,...] 
14 

15 [10,2,1,1,1] [9,4,2] [8,6,1] [5,2,1,...][3,3,2,2,1,...][3,2,2,2,2,2,1,1] 

16 [12,2,1,1][11,5] [11,1,...][8,7,1] [6,1,...][3,2,2,2,2,2,2,1] [4,2,1,...][2,2,2,2,2,1,. 

17 [13,1,...] - - [5,1,...] 

18 - [12,2,1,...][11,4,2,1][10,6,2] [6,2,1,...][4,3,2,2,1,...][3,3,2,2,2,2,1,...] 
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FIG. 1. The density profiles of droplets at v = 1/2 for N = 11 and 17. Filled circles represent the Monte Carlo cal- 
culation of trial wave functions, while the solid lines represent the expected density profile for the case that the angular 
momentum states are evenly occupied (n m = 1/2) up to m e d = 2iV — 1. Note that, although there are some deviations 
between the two, the average value of the Monte Carlo results agree with the expected value. Rwen for the Laughlin states at 
(v,N) — (1/3, 25)(l/5, 20)(l/7, 15) the density profiles have deviations from the expected valueEj. 
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m m 

FIG. 2. Comparison of the result of the exact diagonalization with that of the trial wave function at v = 1/3 for N = 5, 6, 7, 
and 8. Filled circles represent the angular momentum occupation numbers obtained by the exact diagonalization, while the 
open symbols represent these quantities evaluated by a Monte Carlo method using the trial wave functions. 
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m m 

FIG. 3. Comparison of the result of the exact diagonalization with that of the trial wavefunction at v = 1/2 for N = 5, 6, 7, 
and 8. Filled circles represent the angular momentum occupation numbers obtained by the exact diagonalization, while the 
open symbols represent these quantities evaluated by a Monte Carlo method using the trial wave functions. Note that for 
N — 8 there are two possible choices of the trial wave functions, among which [4,2,1,1] has lower Monte Carlo energy. 
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m m 

FIG. 4. The angular momentum occupation numbers for N = 10, 11, 12 and N = 13. Filled circles represent the angular 
momentum occupation numbers obtained by the exact diagonalization, while the open symbols represent these quantities 
evaluated by a Monte Carlo method using the trial wave functions. Among the possible trial wave functions, the lowest energy 
states give consistent results with those of the exact states. 
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FIG. 5. A plot of zeros(o) for a random distribution of electrons(«) for v = 1/2 and N = 12. 
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FIG. 6. The edge occupation numbers are fitted with those of a field theoretical method, Ax + B (solid lines). Filled circles 
and squares represent the exact diagonalization results for N — 11 and N = 12, respectively, and open diamonds represent the 
Monte Carlo result for N = 13 of the trial wave function [5,3,3,1,1] 
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